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ABSTRACT 

As an alternative explanation of the entropy excess and the steepening of 
the X-ray luminosity-temperature relation in groups and clusters, the galaxy 
formation-regulated gas evolution (GG) model proposed recently by Bryan 
makes an attempt to incorporate the formation of galaxies into the evolution 
of gas without additional heating by nongravitational processes. This seems to 
provide a unified scheme for our understanding of the structures and evolution 
of both galaxies and gas in groups and clusters. In this paper, we present an 
extensive comparison of the X-ray properties of groups and clusters predicted by 
the GG model and those revealed by current X-ray observations, using various 
large data sources in the literature and also taking the observational selection 
effects into account. These include an independent check of the fundamental 
working hypothesis of the GG model, i.e., galaxy formation was less efficient in 
rich clusters than in groups, a new test of the radial gas distributions revealed 
by both the gas mass fraction and the X-ray surface brightness profiles, and an 
reexamination of the X-ray luminosity-temperature and entropy-temperature 
relations. In particular, it shows that the overall X-ray surface brightness 
profiles predicted by the GG model are very similar in shape, insensitive to the 
X-ray temperature, and the shallower X-ray surface brightness profiles seen at 
low-temperature systems may arise from the current observational selection 
effect. This can be used as the simplest approach to distinguishing between 
the GG model and the preheating scenario. The latter yields an intrinsically 
shallower gas distribution in groups than in rich clusters. 

Subject headings: cosmology: theory — galaxies: clusters: general — X-ray: 
general 



1. Introduction 


It has been realized for a decade that the ratio of gas mass to stellar mass in groups and 
clusters of galaxies is positively correlated with the X-ray gas temperature. This suggests 
a simple scenario in which galaxy formation was more efficient in groups than in clusters 
(David et al. 1990). As a consequence, the stellar (gas) mass fraction should decrease 
(increase) from groups to clusters, resulting in the increasing mass-to-light ratio and X-ray 
luminosity with X-ray temperature. It seems that these predictions have been well justified 
by the optical measurements of the mass-to-light ratios (e.g. Girardi et al. 2000) and 
the X-ray determinations of the gas fractions in clusters (e.g. Mohr, Mathiesen & Evrard 
1999; MME), and in particular by the discovery of the steepening of the X-ray luminosity 
(L x )-temperature (T) relation for groups and clusters (Henry & Arnaud 1991; Edge & 
Stewart 1991; David et al. 1993; Wu, Xue & Fang 1999; Xue & Wu 2000 and references 
therein). Nevertheless, despite these apparent successes, such a simple scenario has received 
less attention over the past years than the prevailing ‘preheating’ hypothesis. The latter 
assumes that the gas was preheated by nongravitational processes such as supernovae and 
AGNs before collapsing into groups and clusters, which raises the gas entropy and makes 
the warm/hot gas harder to compress (Kaiser 1991). Indeed, the nongravitational heating 
mechanism can significantly reduce the gas content bound within the virial radii of groups 
and poor clusters, providing a satisfactory explanation of the observational facts listed 
above (e.g. Cavaliere, Menci & Tozzi 1997; 1999; Bower et al. 2001). In particular, the 
subsequent detection of the entropy excess in the centers of groups and clusters by Ponman, 
Cannon & Navarro (1999) seems to give a convincing support to the preheating hypothesis. 
Yet, the main difficulty with the preheating model is the unreasonably high efficiency 
of energy injection into the intragroup/intracluster gas from supernovae (Wu, Fabian & 
Nulsen 2000; Bower et al. 2001), although additional heating from AGNs might help to 
alleviate the discrepancy. Furthermore, a uniform preheating of the cosmic baryons to a 
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temperature of ~ 10 6 at high redshifts may make the Lya forest (T ~ 10 4 K) to disappear. 

By inheriting the observational evidence that the efficiency of galaxy formation was 
higher in groups than in clusters, Bryan (2000) proposed that the lowest entropy gas in 
groups and clusters should be responsible for the formation of stars. Namely, the gas in the 
central regions of groups and clusters was converted into stars in the early phase of structure 
formation, and the central gas cavity was refilled with the remaining gas distributed 
originally at large radii. This latter process brought the higher entropy gas inward and thus 
raised the central entropy of the gas, providing an alternative explanation of the existence 
of the entropy floor in the centers of groups and clusters reported by Ponman et al. (1999). 
The most important issue behind this simple scenario, we believe, is that the distribution 
and evolution of the gas are incorporated with the formation of the galaxies so that the 
optical/X-ray properties of the gas and galaxies in groups and clusters can be naturally 
explained within a unified scheme and without invoking additional heating processes. This 
galaxy formation-regulated gas evolution (GG) model in groups and clusters has several 
unambiguous predictions that can be easily tested with current available data or future 
observations. For example, (l)while the stellar and gas mass fractions vary from groups to 
clusters, the total baryon (stellar mass + gas mass) fraction in any group or cluster remains 
a universal value; (2)The mass-to-light ratio is an increasing function of X-ray temperature; 
(3)Less massive systems (e.g. groups) contain less gas and hence have lower X-ray emission, 
which is the key issue for the reproductions of the observed L x -T relation and entropy 
distribution; (4)Radial profile of the stellar mass fraction should demonstrate a decreasing 
trend towards large radii, and an opposite situation occurs for that of the gas mass fraction. 
Of course, whether or not these properties are in quantitative agreement with observations 
deserves further investigation. 

In the pioneering work of Bryan (2000), an overall agreement between the predicted 
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and observed L x -T relations and entropy distributions has been essentially found, 
although at low temperature T < 1 keV, the observed L x -T relation falls below the 
theoretical prediction. Motivated by these successes and their profound implications for our 
understanding of the formation and evolution of the galaxies and gas in groups and clusters, 
we wish to conduct an extensive comparison of the theoretically predicted properties of 
the intragroup/intracluster gas by the GG model with those revealed by current X-ray 
observations. These include not only the global dependence of the X-ray luminosity, entropy 
and gas fraction upon the X-ray temperature T, but also the internal structures of the 
intragroup/intracluster gas such as the radial distribution of gas fraction, the X-ray surface 
brightness profile (slope and core radius parameters), etc. Moreover, some selection effects 
of current X-ray observations (e.g. the finite extensions of X-ray surface brightness profiles 
due to background noise, the spectral fitting temperature, etc.) will be taken into account. 
It is hoped that our comparison would be useful for an eventually conclusive test of the GG 
model. Throughout this paper we assume a flat cosmological model (ACDM) of 12 m = 0.3, 
Ha = 0.7 and h = 0.65. 

2. The Model 

2.1. Reexamination of the working hypothesis 

We begin with a reexamination of the fundamental relationship between the stellar 
mass fraction / star and X-ray temperature T of groups and clusters. Based on three 
independent studies of 33 groups and clusters in the literature, Bryan (2000) derived the 
following empirical formula, regardless of the large scatter of the data points especially for 
low-temperature groups, 


/star = (0.042 ± 0.003) 


kT 

lOkeV 


— 0 . 35 ± 0.06 


( 1 ) 
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where the quoted errors are 68% confidence limits. The lack of a large sample of 
optical/X-ray groups and clusters in which both stellar mass and gas mass are reliably 
determined out to virial radii r vlT hinders us from directly reexamining the validity of 
this relation. Nevertheless, there exist several indirect ways that may allow us to test the 
consistency between equation (1) and other independent measurements. 

The GG model admits the universality of the baryon fraction ft, in different groups and 
clusters: % = / star + / gas . For the ACDM model and = 0.020h -2 , we have j\, = 0.16. 
Therefore, the measurement of the gas fraction / gas in different groups and clusters permits 
an indirect test of the dependence of the stellar mass fraction upon the X-ray temperature. 
The major difficulty of this exercise is the determination of total group/cluster mass M v j r , 
in which one often needs to extrapolate the available X-ray data within finite detection 
aperture to virial radii r vir , along with the assumption of hydrostatic equilibrium. Such an 
extrapolation may be safely applicable only to X-ray luminous and rich clusters. For this 
reason, we extract a sample of 45 nearby (z < 0.1), X-ray luminous clusters from MME 
and Ettori & Fabian (1999; EF). The restriction on cluster redshifts is consistent with 
the three samples used by Bryan (2000), and also gets rid of the influence of cosmological 
parameters other than the Hubble constant. These authors have explicitly derived the gas 
fraction within r 500 where the overdensity A of the cluster dark matter with respect to the 
critical density of the universe is A = 500. rsoo can be approximately regarded as the virial 
radius if we notice that the asymptotic gas fractions of rich clusters at large radii remain 
roughly constant for the conventional isothermal f3 model with /3 = 2/3, although the 
virial radius of a dynamically-relaxed system in the ACDM model corresponds to A & 100. 
We display in Figure 1 the observationally determined gas fraction and the expectation 
°f /gas = 0.16 — 0.042(/cT/10 keV) -0 ' 35 from the empirical formula equation (1) by Bryan 
(2000). It is apparent that the two results are essentially consistent with each other. 



EDITOR: PLACE FIGURE 0 HERE. 


Another independent check of equation (1) is to examine the dependence of the 
mass-to-light ratio M/L upon the X-ray temperature. It follows that the mass-to-light 
ratio M/L can be rewritten as M/L = (M/M star )(M star /L) = T/ s 7 ar , where M star is the 
total stellar mass of groups/clusters and T = M star /L, while observations suggest that 
T « 6.5T 0 for the dominant population of galaxies in clusters, ellipticals (see Fukugita, 
Hogan & Peebles 1998). As a result, the mass-to-light ratio should increase with the X-ray 
temperature according to equation (1). Now, instead of examining the dependence of M/L 
upon T, we would like to derive the f s tax~T relation using the best-fit relation between 
cluster virial mass and optical luminosity found by Girardi et al. (2000): 



( 2 ) 


where c = —1.41 and d = 1.32 for the fore/background correction based on mean galaxy 
counts. This yields 



( 3 ) 


Employing the cosmic virial theorem at £ = 0 (e.g. Bryan & Norman 1998) 



we have 


/star = (l.93 X 10 16 ) 1+1/d lQ-o/d A (l-l/d)/2 



kT 


3(-l+l/d)/2 


lOkeV 


( 5 ) 


Finally, 



( 6 ) 


in which we have adopted T = 6.5T 0 and A = 200. Actually, one can arrive at the 
essentially same conclusion using the empirical M V j r -T relationship for clusters instead of 












the cosmic virial theorem equation (4). For example, the best-fit M 2 oo~T relation based on 
the spatially resolved temperature profiles in the estimates of X-ray cluster masses given by 
Horner, Mushotzky & Scharf (1999) is 


Mo on = 8.2 x 10 14 M P 


kT 


1.48 


0 


lOkeV 


This yields, when combined with equation (3), 


/star = 0.028 


kT 


-0.36 


lOkeV 


( 7 ) 

( 8 ) 


Except the slightly lower amplitudes, both equations (6) and (8) are consistent with 
equation (1) especially when the large uncertainties in the fitting results [equations (1), (2) 
and (7)] are considered. In Figure 1 we illustrate the resulting gas fraction / gas = 0.16 — / star 
using the stellar mass fraction of equation (6). It turns out that the derived gas fraction of 
rich clusters from / star roughly agrees with that determined from the X-ray observations. 


2.2. Basic equations 


The dark matter distributions in groups and clusters are assumed to be unaffected by 
the formation of galaxies and follow the universal density profile (Navarro, Frenk & White 
1995; NFW): 


Pdm(V) = 


^critPcrit 


(9) 


(r/r s )( 1 + r/r s ) 2 ’ 

where 4 cr i t and r s are the characteristic density and length, respectively, and p cr i t is the 
critical density of the universe. The virial mass is defined as M vir = 4wr7 r A/9 cr ; t /3, and 
related to the virial temperature T vir through the cosmic virial theorem equation (4). 
Following Bryan (2000), we adopt the fitting formula of the M vir -c relation suggested by 
numerical simulations: c = 8.5(M v i r h/10 15 M 0 )~ a086 , where c = r v i r /r s is the concentration 
parameter. Since the temperature in all the empirical formulae including the / S tar-F 



9 


relation is the so-called spectral temperature T s rather than the virial temperature T vir , 
an appropriate connection between T s and T vir should be established in order to make a 
meaningful comparison to observations. Using hydrodynamic cluster simulations, Mathiesen 
& Evrard (2001) have recently conducted an extensive comparison of different measures of 
the intracluster gas temperature including the spectral temperature, the emission-weighted 
temperature and the mass-weighted (virial) temperature. In particular, it has been shown 
that T s is generally lower than T vir for clusters, depending on the bandpass used for spectral 
fitting. We adopt the best-fit relation between T s and T v j r in the 2.0-9.5 keV band to 
proceed our numerical computation: T vir = 0.91T S 110 . We have also tested another fitting 
formula in the 0.5-9.5 keV band, T vir = 1.17T S 100 , and found that our main theoretical 
predictions remain unchanged. We will omit the subscript ‘s’ in the spectral temperature. 
Whenever our theoretical predictions are compared with X-ray observations, we always 
work with the spectral fitting temperature. 

For simplicity, we assume that the gas is dissipationless before galaxy formation so 
that the gas traces the underlying dark matter distribution and satisfies the equation of 
hydrostatic equilibrium 

dP° 0 GM DM (r ) 

~T = Pgas v ) - 2 - ; 

ar r 

PgasO) = fbpDM(r). (11) 



We use the superscript ‘0’ to denote the quantities before the formation of galaxies in the 
system. The temperature profile can be obtained by straightforwardly solving the above 
equations with the boundary restriction T°(r —> oo) —» 0: 


T / T 

kT°(r) = kT*— (1 + - 
r* V r s 


(1 + x) ln(l + x) — x 
x 3 (l + x) 3 


where kT* = AnGnm p 5 c p cr i t r^ and /i = 0.59 is the mean molecular weight. Note that unlike 
Bryan (2000), we did not choose the pressure-free external boundary condition P° = 0, i.e. 
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T {] = 0 at r = r vir , which would lead the entropy distribution S° = In T°/{n° e ) 2 ^ to turn 
over near the boundary. While our boundary restriction T°(r —> oo) —> 0 seems to be more 
natural, the validity of equation (10) beyond r vir becomes questionable. In this regard, our 
requirement of T°(r —» oo) —> 0 should not be taken too literally, and the key point here is 
that T°(r) should be a decreasing function of radius near r v i r and tend towards zero when 
r » r vir . 


Galaxy formation is equivalent to removing the central region of radius r 0 an amount 
of gas to M star = f staT M vir : 

/ i-t \ ~ 0 - 35 

M»Jr 0 ) = 0.042 ( IFj - v ) M vlr . 

The remaining gas is redistributed by the conservation of gas mass 



/r 0 


rr 

47 Tn e m p n° e {x)x 2 dx = / 47r /j, e m p n e (x)x 2 dx, 

J 0 


(14) 


where fi e = 2/(1 + A"), and A" = 0.768 is the hydrogen mass fraction in the primordial 
abundances of hydrogen and helium. Namely, the central region of radius r is refilled 
with the gas which is originally distributed at larger radii between r 0 and f. This process 
transports the gas at f with entropy S°(f) = In T°(f)/[n^(f)] 2 ^ 3 to a new position r by the 
conservation of entropy 

T(r) T°(f) 


T°(r 


[n e (r)] 2 / 3 [?7,g (r) ] 2 / 3 ' 

The original position f can be thus fixed by combining the above two equations. Eventually, 
the saturated configuration should satisfy the simple argument 


M gas (r yiT ) + M star (r vir ) = f b M vir . (16) 

Once the equation of state (equation [15]) and the saturated configuration (equation [16]) 
are specified, we can obtain the newly established gas distribution by solving the following 
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equations: 


dP e 


dr 



( 17 ) 


(18) 


in which P e = n e kT. The resulting gas density, temperature and entropy distributions 
are shown in Figure 2 - Figure 4 for various temperatures ranging from 0.5 keV to 14 
keV. It appears that both the cores in the gas density profiles and the entropy floors are 
created in the centers as a result of galaxy formation which has consumed the central gas 
by converting it into stars. The effect is more significant in lower temperature systems than 
in higher ones, as was naturally expected. On the other hand, the conservation of entropy 
during the formation of galaxies leads to a remarkable increase of the gas temperature 
towards the centers of groups/clusters. Although radiative cooling has not been included 
in the above treatment, the significantly raised gas temperature in the central regions can 
help to increase the cooling time which scales as T l P. It thus remains to be an interesting 
issue of whether the combined effect can resolve or release the cooling flow crisis claimed by 
a number of recent observations with Chandra and XMM (e.g. Fabian et al. 2001). For 
the latter, the most compelling evidence comes from the X-ray spectra obtained towards 
the central regions in clusters which show a remarkable lack of emission lines from gas 
with T < 1 keV (Peterson et al. 2001). One of the suggested mechanisms is that there is 
an additional transport mechanism that can reheat the cool gas back to the hot phase. In 
this regard, the formation of galaxies in the central regions which simultaneously raises the 
temperature of the gas may serve as a natural heating mechanism. Of course, further work 
should be done to test this speculation. 


EDITOR: PLACE FIGURE @ HERE. 
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EDITOR: PLACE FIGURE | HERE. 

3. Comparison to observations 
3.1. Radial distribution of gas fraction 

Within the framework of the GG model a significant amount of gas in the central 
regions was converted into galaxies. This would inevitably lead to an increasing gas 
fraction / gas with group/cluster radius, which should be directly testable by current X-ray 
observations. For this purpose, we first extract a sample of 163 clusters, which contain 176 
data points of / gas measured at different radii, from an X-ray image deprojection analysis 
made by White, Jones & Forman (1997). Note that for some of the clusters the X-ray 
temperatures T were estimated by the empirical velocity dispersion-temperature correlation. 
The data points are now divided into three subsamples according to temperature, which 
contain 42 (T > 6 keV, (T) = 7.99 keV), 98 (3 keV < T < 6 keV, (T) = 4.32 keV) and 
36 (T < 3 keV, (T) = 2.36 keV) clusters, respectively, and each subsample is properly 
binned for illustrative purpose. In Figure 5 the radial variations of the observed / gas are 
compared with the predictions for three choices of T = 8, 4.3 and 2.4 keV. Next, we choose 
the two group samples of Mulchaey et at. (1996) and Hwang et al. (1999) used in the 
construction of the / S tar-U relation (see equation [1]) to demonstrate the radial variation of 
the gas fraction in low-temperature systems. These two samples contain 21 groups with 
temperature ranging from 0.69 keV to 3.38 keV, and the mean temperature is (T) = 1.32 
keV. The binned data are displayed in Figure 5. It is immediately clear from Figure 5 that 
there is a fairly good agreement between the model predictions and the observations. The 
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lowest temperature subsample ((T) = 2.36 keV) of White et al. (1997) is an exception, 
exhibiting a somewhat lower amplitude than the theoretical expectation. Perhaps, this can 
be regarded as the most natural explanation of the monotonically increasing gas fraction 
with radius seen in many X-ray clusters although a detailed analysis should be made for 
each individual case. 


EDITOR: PLACE FIGURE | HERE. 


3.2. X-ray surface brightness profile: slope and core radius 


The X-ray surface brightness profile S x (r) can be computed straightforwardly in terms 


of 


SJr) = 


e(T,n e )- 


RdR 


(19) 


2-7T (1 + z) 4 Jr 

where z is the cluster redshift, and e(T,n e ) is the emissivity which is computed by the 
Raymond & Smith (1977) model with a metallicity of Z — 0.3 Z Q . The resulting S x (r) in 
the 0.5-2.0 keV band for nearby clusters (z = 0) are plotted in Figure 6 for a set of 11 
groups/clusters with temperatures ranging from 0.5 to 14 keV. It appears that the overall 
X-ray surface brightness profiles resemble the conventional f3 model in shape. In order to 
facilitate a comparison with the existing X-ray imaging measurements of groups/clusters, 
we fit the predicted X-ray surface brightness profiles to the f3 model characterized by 
S x (r)/S x ( 0) = (1 + r 2 /r;?)U 3/3+1 / 2 ), and work out the slope and core radius parameters, /3 
and r c . Since the determinations of (3 and r c depend somewhat on the extension of the 
X-ray surface brightness profile if the the available data extend only out to two or three 
times the core radius, we truncate each S x (r) at the maximum radii set by the X-ray surface 
brightness limits Sn mit = 2 x 1CT 14 erg s^ 1 arcmin -2 cm -2 and = 2 x 10 -15 erg s -1 
arcmin -2 cm -2 in the 0.5-2.0 keV band, respectively. These two values have approximately 
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covered the ROSAT limits used in current observations of clusters (e.g. MME). It turns out 
that all the predicted X-ray surface brightness profiles can be nicely fitted by the f3 model, 
and the best-fit (3 and r c values for the two flux limits and different choices of redshifts are 
plotted against the X-ray temperature in Figure 7 and Figure 8, respectively. 

EDITOR: PLACE FIGURE | HERE. 

EDITOR: PLACE FIGURE 0 HERE. 

EDITOR: PLACE FIGURE | HERE. 

For a comparison with real observations, we extract a sample of 93 clusters and 36 
groups from the literature in which the best-fit (3 and r c values are explicitly provided. The 
major sources of references include: for clusters, Rizza et al. (1998), EF, MME, Arnaud & 
Evrard (1999), Neumannn & Arnaud (1999), Vikhlinin, Forman & Jones (1999) and Jones 
& Forman (1999); for groups, Mulchaey et al. (1996), Dahlem & Thicring (2000), Helsdon 
& Ponman (2000) and Takahashi et al. (2000). In most cases, the fittings were performed 
by excluding the central cooling components, which is compatible with our model without 
the inclusion of the radiative cooling. The corresponding spectral temperature data are 
taken from the compilation of Wu et al. (1999) and Xue & Wu (2000). 

The theoretically predicted and observationally determined slope and core radius 
parameters for different X-ray temperatures are compared in Figure 7 and Figure 8, 
respectively. Regardless of the large dispersion in the observationally determined (3 and r c 
especially for groups, the trend of the slowly increasing (3 and r c with T is clearly seen, 
which is consistent with our predictions. In particular, the predicted slope parameter varies 
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from [3 ~ 0.3 for groups to (3 ~ 0.9 for very rich clusters, in excellent agreement with current 
observations. It should be pointed out that the trend towards shallowed surface brightness 
profiles in low-temperature systems can be attributed to an observational selection effect 
rather than a property intrinsic to groups/clusters. This can be clearly seen from the 
overall predicted X-ray surface brightness profiles for systems with different virial masses. 
Such a property differs significantly from the preheating model in which the flattening of 
the S'x occurs in a low-temperature system because the preheated gas becomes harder to be 
trapped in the shallow gravitational potential well (see Cavaliere et al. 1997). We have also 
plotted in Figure 7 the slope parameters expected in two extreme situations for comparison 
(cf. Appendix): (l)The isothermal model in which the gas is driven by purely gravitational 
shocks and (2) the isentropic model in which the gas is preheated and then collapses 
adiabatically into groups/clusters. The former may correspond to rich clusters, while the 
latter may be applied to groups under the preheating hypothesis. These two extreme cases 
give approximately rise to the upper and lower bounds to the (3 values. 

Whether the increasing slope parameters [3 of X-ray surface brightness profiles from 
groups to cluster seen in current observations is purely an observational selection effect or a 
property intrinsic to groups/clusters constitutes a crucial test for the GG model. There are 
at least two major uncertainties in the current f3 model fits of the X-ray surface brightness 
profiles of groups: Firstly, for most groups a single f3 model does not provide an adequate 
description of the observed data, and the fit quality depends sensitively on the modeling of 
central galaxies, cooling flows or AGNs (e.g. Helsdon & Ponman 2000); Secondly, the X-ray 
emission of most groups remains undetectable outside radii of typically ~ 200-400 kpc clue 
to the limitation of current detector sensitivity (e.g. Ponman et al. 1996; Mulchaey et al. 
1996; Helsdon & Ponman 2000), while the virial radii of galaxy groups with temperature 
T 1 keV are about 1 Mpc. It seems that current X-ray observations have only probed 
the central regions other than the whole groups. It is hoped that the high sensitivity 
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measurements with the Chandra X-ray Observatory can provide a robust constraint on the 
slopes of X-ray surface brightness profiles for groups and poor clusters especially at large 
radii. 


3.3. X-ray luminosity - temperature relation 

We now come to the global properties of groups and clusters. The first test that Bryan 
(2000) chose for the GG model is the X-ray luminosity-temperature relation. While the 
predicted L x -T relation is essentially consistent with the observed data on cluster scales, 
there is an apparent disagreement below T « 1 keV. This could be due to the small sample 
of poor clusters and groups and/or the observational selection effect of the finite spatial 
extensions of X-ray emission set by the X-ray flux-limited surveys. For the latter possibility, 
we recall that the correction for lost flux falling out the detection aperture has been already 
made in the computation of L x for clusters. The situation for the groups with (3 < 1/2 is 
complicated. In most cases, no aperture correction is made in order to avoid the divergence 
of total X-ray luminosity (e.g. Ponman et al. 1996; Helsdon & Ponman 2000). It is thus 
possible that the X-ray luminosity excess on group scales predicted by Bryan (2000) can 
be reduced simply by excluding the X-ray emission outside the detection aperture implied 
by the flux limit We now re-examine the L x -T relation using the catalog of X-ray 

groups and clusters compiled by Wu et al. (1999) and Xue & Wu (2000). The updated 
catalog contains 57 groups and 192 clusters whose X-ray temperatures and luminosities are 
both available. Furthermore, we only account for the X-ray emission within the apertures 
set by the X-ray surface brightness limits »S'i imit = 2 x 10” 14 erg s -1 arernin -2 cm -2 and 
S'limit = 2 x 10- 15 erg s' 1 aremin -2 cm -2 in the 0.5-2.0 keV band, respectively. It appears 
that the evaluation of the total X-ray luminosity of clusters is almost unaffected by these 
aperture limits because the corresponding cutoff radii are close to or beyond the virial 
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radii. For the latter, we simply truncate the clusters at their virial radii. Again, we use the 
Raymond & Smith (1977) model with a metallicity of Z = 0.3 Z Q to calculate the X-ray 
luminosity and plot the resulting L x -T relation in Figure 9. It turns out that the model 
matches nicely the observed data over entire temperature range. 

EDITOR: PLACE FIGURE | HERE. 


3.4. Entropy distribution 

The last check is the distribution of the central entropy S measured at 0.1r 2 oo against 
the X-ray temperature. We use the updated measurements of S'(O.lrvir) by Ponman et al. 
(1999), Lloyd-Davies, Ponman & Cannon (2000) and Xu, Jin & Wu (2001), which contain 
67 data points ranging from T fa 0.5 keV to T fa 14 keV. We demonstrate the predicted 
and measured S^O.lrVir) in Figure 10, which reinforces the finding of Bryan (2000) that 
the predicted S-T is consistent with the observed one. A larger sample especially at low 
temperature end will be needed for a critical test of the model prediction. 


EDITOR: PLACE FIGURE ITTl HERE. 


4. Discussion and conclusions 

Using various large data sources of X-ray groups and clusters in the literature, we have 
made an extensive examination of the properties of the galaxy formation-regulated gas 
evolution model in groups and clusters, proposed recently by Bryan (2000) as an alternative 
explanation of the entropy excess and the steepening of the X-ray luminosity-temperature 
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relation in groups and clusters. This model is based on one empirical correlation between 
the stellar mass fraction / star and temperature T and one hypothesis that the lowest entropy 
gas in the central regions of groups and clusters was converted into stars. The reliability of 
the /star-7" relation has been justified by other observational facts such as the increasing 
ratio of gas mass to stellar mass with the X-ray temperature (David et al. 1990), the 
increasing gas fraction with the X-ray temperature (Figure 1) and the positive correlation 
between the mass-to-light ratio and the optical luminosity in clusters (Girardi et al. 2000). 
This suggests a simple scenario that galaxy formation was less efficient in rich clusters than 
in poor clusters and groups. As for the working hypothesis, the formation of galaxies in the 
central regions of groups and clusters consumed the central gas and created a gas cavity, 
which acts as a trigger for an inward gas flow or diffusion. As a consequence, the flat gas 
core and entropy floor developed in the centers, together with an increasing gas fraction 
with radius. These predictions are in quantitative agreement with the internal structures 
of intragroup/intracluster gas revealed by current X-ray observations such as the slope 
and core distributions characterized by the conventional (5 model, the radial variations of 
the gas mass fraction, and the central entropy distributions measured at 0.1r 2 oo- Another 
success of the model is the reproduction of the X-ray luminosity-temperature relation over 
a broad temperature range from T « 0.5 keV to T & 20 keV. Overall, it seems that this 
simple model without additional heating can satisfactorily account for the existing X-ray 
observations. 

While both the GG model proposed by Bryan (2000) without additional heating and 
the prevailing preheating model can equally explain the current X-ray observations, there 
are a number of essential differences which can be used to distinguish between the two 
scenarios. For example, the GG model admits the universality of the baryon fraction from 
groups to rich clusters, while in the preheating scenario the baryon fraction increases with 
temperature because a considerable amount of gas may still reside outside of the virial radii 
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as a result of preheating especially in poor clusters and groups. A precise measurement of 
the total baryon (stars + gas) fractions in groups and clusters can provide a critical test for 
the models. Alternatively, the X-ray imaging measurements of the surface brightness profiles 
of poor clusters and groups are probably the simplest way to disentangle the issue: the GG 
model predicts that the slope of the X-ray surface brightness profile at large radii remains 
roughly unchanged from groups to clusters, and the increasing slope with temperature seen 
in current X-ray observations is purely selection effect due to the limitation of detector 
sensitivity and background noise. On the contrary, within the framework of the preheating 
scenario the X-ray surface brightness profiles are intrinsically shallower in groups and poor 
clusters than in rich clusters, insensitive to observational sensitivity. The high sensitivity 
X-ray observations of groups with advanced detectors like Chandra should be able to set 
robust constraints on the slope parameters of the X-ray surface brightness profiles in the 
outer regions of groups. 

Apparently, the GG model still has its own problems that need to be resolved in the 
future. First, aside from the robustness problem of the empirical relation between stellar 
mass fraction and X-ray temperature established based on a small sample of local groups 
and clusters, it is unclear as to how this relation evolves with redshift. The cosmological 
applications are largely limited if evolutionary effect is not incorporated into the present 
model. It is unlikely that there is a significant evolution of groups and clusters since z ~ 1 
if the inward gas flow or diffusion due to galaxy formation in the central regions ended at 
the early phase of structure formation. Of course, this depends on how fast the process of 
inward gas flow or diffusion may take. A related question is: Can one actually observe the 
process at high-redshift groups and clusters ? Second, the current / sta r-T relation indicates 
that nearly all the gaseous materials in less massive groups and giant galaxies have been 
converted into stars. If this occurred in the early stage of structure formation, there would 
be no gas left in small systems. How can one explain the existence of a considerably large 
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fraction of hot gas in clusters if they formed by gravitational aggregation of individual 
low-mass objects ? One possibility is that the formation of galaxies and the inward gas flow 
or diffusion in less massive systems proceeded rather slowly and a large fraction of gas had 
not been converted into stars before they merged with other massive objects. The latter 
process can raise the gas temperature and thus reduce the efficiency of galaxy formation, 
according to the / star -T relation. A sophisticated analysis incorporated with the halo 
merger rates described by the extended Press-Schechter formalism (Lacey & Cole 1993) 
may be needed to clarify the issue. Third, the GG model predicts that the temperature 
profile shows a dramatic increase towards the centers and a significant decline at large radii. 
Although the inclusion of radiative cooling is expected to reduce somewhat the central 
X-ray temperature, whether or not this can yield a quantitative agreement with X-ray 
spectral analysis remains unclear. It is also interesting to explore the possibility of whether 
or not such a rising temperature towards the centers of clusters can alleviate the cooling 
flow crisis. Finally, it is worth mentioning that the overall radiative cooling may play an 
important role in the removal of low-entropy gas in the centers of groups and clusters, 
giving rise to an equally good explanation of the steepening of the L x -T relation for groups 
and clusters (Muanwong et al. 2001). 

Overall, the GG model has seemingly proved very successful at explaining a number of 
internal and global properties of the intragroup/intracluster gas. If confirmed, this would 
have profound implications for our understanding of the origin and distribution of the 
intragroup/intracluster gas, the formation history of group/cluster galaxies, the missing 
baryons in the universe, and even the formation and evolution of structures in the universe. 
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In this Appendix we discuss the gas distributions in two extreme cases: the isothermal 
and isentropic models. The former may work in rich clusters where the gas was mainly 
heated by gravitational shocks, while the latter may correspond to the intragroup gas which 
was preheated and then collapsed adiabatically into the gravitational potential wells of the 
groups. A detailed discussion about the X-ray luminosity-temperature relations in the two 
models can be found in Balogh et al. (1999). Here, we concentrate on the radial gas density 
profiles and use a new prescription suggested recently by Eke, Navarro & Steinmetz (2001) 
for the assignment of the collapse redshift of dark halos. Again, we adopt the NFW profile 
for a virialized dark halo, and determine its virial temperature in terms of equation (4). In 
order to fix the free parameters, 8 C or r s or c in the NFW profile. The collapse redshift z co \\ 
for each halo identified at z = 0 is introduced through (Eke et al. 2001) 

B(z C oii)cr eff (M s ) = (1) 

a 


where C a ~ 25 for the ACDM model, D is the normalized linear grow factor, for which 
we take the approximate expression from Carroll, Press & Turner (1992), a e s(M s ) is the 
so-called modulated rms linear density at mass scale M s : 


Vett(M s ) 


a(M s ) 


d\i\cr(M s ) 
d In M s 


( 2 ) 


and M s is defined as the mass contained within r max = 2.17r s where the circular velocity 
reaches its maximum and 


AL = 


0.47M 


ln(l + c) — c/(l + c)‘ 

We adopt the following parameterization of the power-spectrum of initial fluctuation 


( 3 ) 


P(k) = A k n T* DM (k), 


( 4 ) 


where n is the primordial power spectrum and is assumed to be the Harrison-Zeldovich case 
n — 1, and Xcdm(^) is the transfer function of adiabatic CDM model for which we use the 
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fit given by Bardeen et al. (1986): 

Tcb m( 9) = ln(1 2 + 3 4„ 34<?) l 1 +3.89 9 + (16.1 g 2 ) + (5.46<?) 3 + (6.71< ? ) J ]- 1 / 4 , (5) 

where q = ( k/h Mpc _1 )/r, and T = ^M^exp[—flb(l + V^h/^t m)] is the shape parameter. 
Once a power spectrum, P(k), is specified, the mass variance becomes 

1 r°° 

cr 2 (Ad) = — / k 2 P(k)W 2 (kR)dk, (6) 

Z7 T 2 Jo 

in which W(x) = 3 (sin re — xcosx)/x 3 is the Fourier representation of the window function. 
The amplitude A in the power spectrum is determined using the rms mass fluctuation 
on an 8 h~ l Mpc scale, as, for which we take a§ = 0.93 for our ACDM model. Finally, 
the relationship between the NFW profile and its collapse reshift z co n is established if we 
define a new characteristic density, p s = 3M v i r /47rrf, and set to equal the spherical collapse 
top-hat density at z co u: 

Ps ^(^zoll)Pcrit (^coll) • if) 


I. Isothermal model 


We first assume that the hot gas is isothermal with electron number density n e (r ) and 
temperature T, and is in hydrostatic equilibrium with the underlying gravitational potential 
dominated by the NFW profile: 

_ GAd DM (r) _ 1 d[kTn e (r )] 

r 2 pm p n e (r) dr 

A straightforward computation yields an analytic form of the electron density profile 
(Makino, Sasaki & Suto 1998) 


n e (r ) 


n e ( 0) 


(1 + r/r s )“ /(r/r#) 

e Q - 1 



(9) 


where we have adopted the normalized, background-subtracted form in order to ensure the 
convergence of the X-ray surface brightness, and a = AnGimi p 8 c p cr [tr 2 /kT. In terms of the 
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virial theorem, the a parameter reduces to 

« = —--V (10) 

ln(l + c) - c/(l + c) 

Here we have assumed that the specific kinetic energy of dark matter particles is equal to 
that of gas: /3 spec = cr 2 / (kT vil /pm p ) = 1. In turns out from the prescription of the dark 
halo formation discussed above that a is a slowly decreasing function of virial mass or 
temperature, ranging from 22.5 to 14 for 1O 9 M 0 < M < 1O 17 M 0 . This is in good agreement 
with the observationally determined values for clusters (e.g. EF; Wu & Xue 2000). 

Once the electron density and temperature are specified for different systems, we can 
calculate the X-ray surface brightness profile using an optically thin, isothermal plasma 
emission model by Raymond & Smith (1977). Finally, we fit the predicted X-ray surface 
brightness profile to the f3 model and work out the best-fit /3 and r c parameters. It appears 
that the resulting f3 parameter decreases slightly with temperature, as shown in Figure 
7. At the high temperature end of T > 10 keV, the predicted slope parameter reaches 
f3 ~ 0.8 — 0.9, roughly consistent with observed values for very rich clusters. 


II. Isentropic model 


The specific entropy of the gas can be conveniently defined as 


S 


kT 


7—1 

Tie 


i 


( 11 ) 


where 7 is the polytropic index. In the case of where the gas is accreted adiabatically, 

7 = 5/3. We solve the hydrostatic equilibrium equation by demanding that the entropy 
of the gas is conserved during accretion. We specify such a boundary condition that the 
electron number density and temperature at large radii should approach asymptotically the 
background values n e ,6 and T&, respectively. We caution that this last constraint may break 
down beyond virial radius because of the failure of hydrostatic equilibrium. Therefore, this 
boundary condition should only be taken to be a reasonable approximation. The electron 
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number density and temperature profiles are 

'7 - l\ ln(l + r/r s ) 


n e (r) = n efi 


1 + Qt b 


1 

7-1 


T(r ) = T b 


1 + 0L b 


7 ) r/r s 
(l ~ l\ ln(l + r/r s ) 


7 


r/r s 


( 12 ) 

(13) 


where 


Oib 


4nGimip5 c p cr it j rs _ 4:7rGpm p 6 c p cr i t r, 


m 


Sn ^ 1 


(14) 


Nevertheless, we use the background-subtracted electron density n e (r) = n e (r ) — n e> b and 
temperature T(r) = T(r) — T b to proceed our calculation of the corresponding X-ray 


emission, i.e., 


n e (r) 
Tle,b 


l+Qf^ 


'7 - l\ ln(l + r/r s ) 


7 


1 

7-1 


1 ; 


T(r) 

T h 


= OL b 


r/r s 

'7 —■ l\ ln(l + r/r s ) 


7 


r/r s 


(15) 

(16) 


The key parameter that determines the shape of the X-ray surface brightness profile S x (r) 
is a b . The typical value of a b in terms of the above evaluations of 8 C and r s for a given 
halo M vir or T vir is a b ~ 4:TTGp/rn p 5 c p rr \trl /A;T vir ~ 10 for 1O 9 M 0 < M < 1O 17 M 0 . We have 
numerically obtained the profiles of S x (r) for a number of a b and fitted them to the f3 model 
over a radius range of 0 < r < 100r s . The resulting f3 starts from ~ 0.3 for a b = 10 and 
approaches ~ 0.4 for the extremely large values of a b ~ 10 4 (see Figure 7). Such a broad 
range of a b should cover a broad mass range from groups to rich clusters. It appears that 
the saturated configuration of the isentropically accreted gas in an NFW-like gravitational 
potential well looks quite shallow and can be well represented by the conventional f3 model 
with /3 0.3 — 0.4. 
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Fig. 1.— Gas mass fraction within r 5 oo versus X-ray temperature for nearby rich clusters. 
A total of 55 clusters from Mohr et al. (1999; MME) and Ettori & Fabian (1999; EF) are 
shown. Also plotted are the expectations from the correlations between the stellar mass 
fraction and X-ray temperature given by Bryan (2000) (solid line) and between the virial 
mass and optical luminosity found by Girardi et al. (2000) (dotted line), in combination 
with the universality of the total baryon fraction in clusters predicted by the Big Bang 
Nucleosynthesis (BBN) (dashed line). 

Fig. 2.— Normalized gas density profiles for different X-ray temperatures ranging from 0.5 
keV to 14 keV. The initial profiles and the galaxy formation-regulated gas distributions are 
plotted by dotted and solid lines, respectively. 

Fig. 3.— The same as Figure 2 but for the true temperature profiles. 

Fig. 4.— The same as Figure 2 but for the entropy profiles. 

Fig. 5.— Radial variations of the gas mass fraction. The predicted profiles (dotted lines, 
from top to bottom) correspond to four different choices of temperatures T = 8.0, 4.3, 2.4 
and 1.3 keV. Also shown are the observationally derived 176 data points of clusters (open 
symbols) by White et al. (1997) and 21 data points of groups (filled stars) by Mulchaey et al. 

(1996) and Hwang et al. (1999) within different radii, which are properly binned according 

to temperature and radius. 

Fig. 6.— The predicted X-ray surface brightness profiles in the 0.5-2.0 keV band for a set 
of 11 groups/clusters with temperatures ranging from 0.5 to 14 keV (from bottom to top). 
Dotted lines correspond to the X-ray surface brightness limits ,S'ij m i t = 2 x 10~ 14 erg s^ 1 
arcmin -2 cm -2 (upper) and A ihn it = 2 x 10~ 15 erg s -1 arcrnin -2 cm -2 (lower), respectively. 

Fig. 7.— Dependence of the slope parameter (/3) upon X-ray temperature (T). A total of 
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127 groups and clusters extracted from the literature are shown. Thick and thin solid lines 
are the best-fit f3 values from the GG model with and without the correction for maximum 
detectable radius, among which the three thick solid lines from top to bottom correspond 
to (S'limit, z)= (2 x 10 -14 erg s' 1 arcmin -2 cm -2 , 0), (2 x 10 -15 erg s -1 arcmin -2 cm -2 , 0), 
and (2 x 10 -15 erg s -1 arcmin -2 cm -2 , 0.2), respectively. Dashed line is the best-fit value 
from the isothermal model. Dotted lines represent the results in the isentropic model for 
ranging from 10 to 10 4 , which covers a broad mass range from groups to rich clusters (see 
Appendix). 

Fig. 8.— Comparison of the theoretically predicted core radii and the best-fit values from 
X-ray observations for different temperatures. Virial theorem is used to assign r 2 oo to each 
group/cluster. Thick and thin lines correspond to the f3 model fittings with and without the 
correction for maximum detectable radius. The notations are the same as in Figure 7. 

Fig. 9.— Comparison of the predicted bolometric X-ray luminosity-temperature relation 
and the observed data. A total of 57 groups and 192 clusters with available T and L x in the 
literature are shown. The upper and lower lines correspond to the two flux limits shown in 
Figure 6, 5 7,-n it. = 2 x 10 -14 erg s -1 arcmin -2 cm -2 and S An it. = 2 x 10 -15 erg s -1 arcmin -2 
cm -2 , respectively. 

Fig. 10.— Dependence of the central gas entropy measured at 0.1r 2 oo upon the X-ray 
temperature. The model prediction is plotted by solid line. The observational results from 
Ponman et al. (1999; PCN), Lloyd-Davies et al. (2000; LPC) and Xu et al. (2001; XJW) 
are clearly marked. 
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